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Abstract 

Within a light-cone QCD formahsm based on the Green function technique incorpo- 
rating color transparency and coherence length effects we study nuclear shadowing in 
deep-inelastic scattering at moderately small Bjorken xbj- Calculations performed so far 
were based only on approximations leading to an analytical harmonic oscillatory form 
of the Green function. We present for the first time an exact numerical solution of the 
evolution equation for the Green function using realistic form of the dipole cross section 
and nuclear density function. We compare numerical results for nuclear shadowing with 
previous predictions and discuss differences. 
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1 Introduction 



Nuclear shadowing in deep-inelastic scattering (DIS) off nuclei is intensively studied during 
the last two decades. It can be treated differently depending on the reference frame. In 
the rest frame of the nucleus this phenomenon looks like nuclear shadowing of the hadronic 
fluctuations of the virtual photon and is occurred due to their multiple scattering inside the 
target 0121131111301111313 CD]- In the infinite momentum frame of the nucleus it can 
be interpreted, however, as a result of parton fusion [TT[ [T21 [13 Ull leading to a reduction 
of the parton density at low Bjorken xbj- Although these two physical interpretations are 
complementary, we will work in the rest frame of the nucleus, which is more intuitive and is 
well suited also for the study of the coherence effects jT3]. 

Important phenomenon which controls the dynamics of nuclear shadowing in DIS is effect 
of quantum coherence. It results from destructive interference of the amplitudes for which the 
interaction takes place on different bound nucleons. It can be treated also as the lifetime of the 
qq fluctuation and estimated by relying on the uncertainty principle and Lorentz time dilation 
as, 

- 2t/ 

where u is the photon energy, is photon virtuality and Mgg is the effective mass of the qq 
pair. It is usually called coherence time, but we also will use the term coherence length (CL), 
since light-cone kinematics is assumed, Ic = t^. CL is related to the longitudinal momentum 
transfer q^ = 1/lc- The effect of CL is naturally incorporated in the Green function formalism 
already applied in DIS, Drell-Yan pair production [1312] and vector meson production pT| [TH] 
(see also the next Section). 

The nuclear shadowing in DIS was studied in ^3 using correct quantum mechanical 
treatment based on the Green function formalism. The Green function controls then not only 
the relative transverse motion of the qq pair but also an importance of the higher order multiple 
scatterings in the nucleus. The solution of the evolution equation for the Green function was 
performed so far analytically. This analytical solution requires, however, to implement several 
approximations into a rigorous quantum-mechanical approach like a constant nuclear density 
function (see Eq. (|7T|)) and a specific quadratic form of the dipole cross section (see Eq. ([20))). 
Consequently, obtained in a such way the harmonic oscillator Green function (see Eq. (|22p) was 
used for calculation of nuclear shadowing. However, the following question naturally arises; how 
accurate is the evaluation of the nuclear shadowing in DIS using this Green function ? In order 
to clarify this one should solve the evolution equation for the Green function numerically. 
It does not bring any additional assumptions and does not force us to use supplementary 
approximations, which cause the theoretical uncertainties. Therefore the main goal of this 
paper is to present for the first time the predictions of nuclear shadowing in DIS at moderately 
small XBj based on exact numerical solution of the evolution equation for the Green function. 
In addition, applying an algorithm described in the Appendix A we present also calculations of 
nuclear shadowing within the harmonic oscillator Green function approach using quadratic form 
of the dipole cross section (Eq. ([13)) and a constant nuclear density function (Eq. ^^). We 
check whether they correspond to the results already presented in JTSJ . Finally we analyze and 
discuss the differences between the exact and approximate predictions for nuclear shadowing. 
Advantages of an exact numerical solution of the two-dimensional Schrodinger equation for the 
Green function (see Eq. fll 7|) ) presented in this paper provide a better baseline for the future 
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study of the QCD dynamics not only in DIS off nuclei but also in further processes occurred 
in lepton (proton)-nucleus collisions. 

Calculations of nuclear shadowing presented in the paper ^3] were performed assuming only 
qq fluctuations of the photon and neglecting higher Fock components containing gluons and 
sea quarks. Performing realistic calculations, we include the effects of higher Fock states as the 
energy dependence of the dipole cross section, crgq(r, s)^. We use two realistic parametrizations 
of (Tqqif^ s) (scc the ucxt Scctiou and Eqs. (0) and (0)), where the energy dependence is naturally 
included. However, we will neglect higher Fock states leading to gluon shadowing (GS) 
assuming only low and medium values of the photon energy v as was done also in [T5] . 

The paper is organized as follows. In the next Section we present the light-cone dipole 
phenomenology for nuclear shadowing in DIS together with the Green function formalism. 
The Section El supplemented by Appendix A is devoted to description of an algorithm for 
numerical solution of the evolution equation for the Green function. Numerical results based 
on realistic calculations and a comparison with predictions within harmonic oscillator Green 
function approach are presented in Section |31 Finally, in Section we summarize our main 
results and discuss differences between realistic and approximate [13 HE] calculations of nuclear 
shadowing in DIS. 



2 Light-cone dipole phenomenology for nuclear shadow- 
ing 

The main goal of the light-cone (LC) dipole approach to nuclear shadowing is a possibility 
to include the nuclear form factor in all multiple scattering terms. Derivation of the formula 
for nuclear shadowing can be found in [201 • The study of the difference between the correct 
quantum-mechanical treatment of nuclear shadowing and known approximations is given in 
ITS' assuming only qq Fock components of the photon and neglecting higher Fock components 
containing gluons and sea quarks. The nuclear antishadowing effect was omitted as well because 
was assumed to be beyond the shadowing dynamics. The total photoabsorption cross section 
on a nucleus can be formally represented in the form 

<y^'^{xBj, Q') = A a^*^(a;Bj, Q') - Aa(xB„ Q^) . (2) 
Here the Bjorken variable xsj is given by 

where s is the 7*-nucleon center of mass (cm.) energy squared and is mass of the nucleon. 
cr^*^ {xBj,Q^) in ^ is total photoabsorption cross section on a nucleon 

Q') = / rfV £ da |^,-,(r, a, Q^) ^ agg{f, s) . (4) 

Here crqq{r, s) is the dipole cross section which depends on the qq transverse separation r and 
the cm. energy squared s and "^qqir, a, Q"^) is the LC wave function of the qq Fock component 

^Here r represents the transverse separation of the qq photon fluctuation and s is the center of mass energy 
squared (see the next Section). 
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of the photon which depends also on the photon virtuahty and the relative share a of the 
photon momentum carried by the quark. Note that Bjorken xbj is related with cm. energy 
squared s via Eq. Q. Consequently, hereafter we will write the energy dependence of variables 
in subsequent formulas also via xsj- dependence whenever convenient. 

The first ingredient of the photoabsorption cross section on a nucleon (j^ is the dipole 
cross section ijqq{r,s) representing interaction of a qq dipole of transverse separation r with 
a nucleon [21 . It is a flavor independent universal function of r and energy and allows to 
describe various high energy processes in an uniform way. It is known to vanish quadratically 
CTqq{r, s) oc as r ^ due to color screening (property of color transparency j2ll |22l I^Hj ) and 
cannot be predicted reliably because of poorly known higher order perturbative QCD (pQCD) 
corrections and nonperturbative effects. There are two popular parameterizations of (Tgg(r, s), 
GBW presented in |^ and KST suggested in [121 . Detailed discussion and comparison of 
these two parametrizations can be found for example in [23 El- Therefore, for completeness, 
we present here only the main features of both parametrizations because they are used in the 
realistic calculations of nuclear shadowing in DIS with the results shown in Section |3] 

The GBW model |21] for the dipole cross section provides a very simple parametrization 
which saturates at large qq separations. 



where Roixsj) = 0.395 (xbj/xo)^/^ fm and ctq = 23.03 mb; A = 0.288; Xq = 0.0003. This dipole 
cross section vanishes oc at small dipole sizes as implied by color transparency (CT). It 
describes well the data for DIS at small xbj and medium and large Q^. However, it cannot be 
correct at small since predicts energy-independent hadronic cross sections. Besides, xsj is 
not any more a proper variable at small and should be replaced by energy. This problem 
is removed by the KST parametrization which keeps the form © but contains an explicit 
dependence on energy. 



An explicit energy dependence in the parameter o"o(s) is introduced in a such way that guar- 
antees the reproduction of the correct hadronic cross sections. 



where crtJl{s) = 23.6 (s/so)°'°^^ + 1.432 (s/sq) °''*^mb is the Pomeron and Reggeon parts of the 
np total cross section and Ro{s) = 0.88 {s/ sq)'^/'^ im with A = 0.28 and sq = 1000 GeV^ 
is the energy-dependent radius. In Eq. ([7j) {r'^h)^ = 0.44 fm^ is the mean pion charge radius 
squared. The main advantage of the KST parametrization ([7j) is that it describes well the 
transition down to limit of real photoproduction, = 0. However, the improvement compared 
to GBW model [24j at large separations (small values of Q^) leads to a worse description of 
the short-distance part of the dipole cross section which is responsible for the behavior of the 
proton structure function at large Q^. To satisfy Bjorken scaling, the dipole cross section at 
small dipole sizes r must be a function of the product sr which is not the case for the KST 
parametrization The form of Eq. ^ successfully describes the data for DIS at small xbj 
only up to ~ 10 GeV^ and does a poor job at larger values of Q^. 








(7) 
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Summarizing, the GBW model is suited better at medium and large > 5 ^ 10 GeV^^ 
and at medium small and small xbj ^ 0.01 whereas the KST model prefers low and medium 
values of < 5 -i- 10 GeV^. Therefore, the difference of the realistic calculations for nuclear 
shadowing in DIS using these two models for the dipole cross section in the common kinematic 
region of their applicability can be treated as a measure of theoretical uncertainty. 

The second ingredient of <j'^*^{xBj,Q^) in (jH) is the perturbative distribution amplitude 
("wave function") of the qq Fock component of the photon^ and has the following form for 
transversely (T) and longitudinally (L) polarized photons j2Zl 123 H] , 



where x and x are the spinors of the quark and antiquark, respectively; Zq is the quark charge, 
Nc = 3 is the number of colors. KQ{er) is a modified Bessel function with 

e2 = a (1 - a) + , (9) 
where niq is the quark mass. The operators O^'^ read, 

= a ■ e + i (1 - 2a) {a ■ n) (e ■ Vr-) + (a x e) ■ V,. , (10) 

0^ = 2Qa{l- a){a -n) . (11) 

Here Vr acts on transverse coordinate r\ e is the polarization vector of the photon, n is a unit 
vector parallel to the photon momentum and a is the three vector of the Pauli spin-matrices. 

Matrix element (0)) contains the LC wave function squared, which has the following form 
for T and L polarization, 

9 AT 

*?,(^«,Q^) f = E [m)K,{e,rf + [a' + (1 - af] K,{erf] , (12) 

and 

n,{r,a,Q') f = ^-j^ E Z}Q'a\l - af K,{erf , (13) 

where Ki is the modified Bessel function, 

K,{z) = -j-K,{z). (14) 

Note that in the LC formalism the photon wave function contains also higher Fock states 
\qqG), \qq2G), etc. The effects of higher Fock states are implicitly incorporated into the 
energy (Bjorken xbj- ) dependence of the dipole cross section agq{r,s) as is given in Eq. (jl)). 
Note, that the energy dependence of the dipole cross section is naturally included in realistic 
parametrizations Eqs. (jSj) and ©. 

i.e. at medium small and small values of dipole size, r oc yT7Q2< 0.06 ^0.09 fm. 

^We neglect the nonperturbative effects responsible for the interaction between q and q assuming sufficiently 
large values of in DIS (see below). 
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In Eq. (j2)) the second term, Act, represents the shadowing correction and has the following 
form 

2 r poo rOG pi 

Aa{xBj,Q'^) = - Re cfb dzipAib,zi) dz2 PAib, Z2) da A{zi, Z2, a) , (15) 
2 J J— 00 J zi Jo 

with 

A{zi,Z2,a) = J rfVs ^*-g(r^,a,(5^) agg{r2,s) J d^ri (^^^(r^, ^2; n, ^1) o-g<y(ri,s) ^^^(rl, a, Q^) . 

(16) 

In Eq. (fT3j) pa(&, z) represents the nuclear density function defined at the point with longitudinal 
coordinate z and impact parameter h. 




Figure 1: A cartoon for the shadowing term Act in (0)- Propagation 
of the qq pair through the nucleus is described by the Green function 
Gqq{f*2, Z2]f\, Zi) which results from the summation over different paths 
of the qq pair. 

The shadowing term in Q is illustrated in Fig. ^ At the point zi the initial photon 
diffractively produces the qq pair (7*A^ — > qqN) with transverse separation fi. The qq pair then 
propagates through the nucleus along arbitrary curved trajectories, which are summed over, 
and arrives at the point Z2 with a transverse separation f2. The initial and final separations 
are controlled by the LC wave function of the qq Fock component of the photon ^^^^(r, a, Q^). 
During propagation through the nucleus the qq pair interacts with bound nucleons via dipole 
cross section (Jqq{r, s) which depends on the local transverse separation r. The Green function 
Gqq{f'2, Z2;f\,zi) dcscribcs the propagation of the qq pair from zi to Z2. 

The Green function Gqq{f^, Z2;'fl,zi) satisfies the time-dependent two-dimensional Schrodinger 
equation, 

Gqq{f*2, Z2; f{, Zl) = 

dZ2 

with the boundary condition 

Ggg(r^,Z2;n,Zi)|22=2i =5^(n-^2). (18) 



A 



r2 



2u ail ~ a 



+ Vqq{Z2,r2,a) 



Gqq{r2,Z2]fi,Zi) 



(17) 



6 



In Eq. (|17p the Laplacian acts on the coordinate r and e is given by 0- 

The Green function Gqq{f^, Z2','fi, zi) includes the phase shift between initial and final pho- 
tons which is due to transverse and also longitudinal motion of the quarks. One can see the 
presence of the CL in the kinetic term of the evolution Schrodinger equation (|T7j) . where 
the role of time is played by the longitudinal coordinate Z2. A part Ar/[2z/a(l — a)] of this 
kinetic term takes care of the varying effective mass of the qq pair, M|g = (m^ + k^)/a{l — a), 
and provides a proper phase shift. This is what the overall kinetic term consists when the 
transverse momentum squared of the quark is replaced by — > A^. This dynamically varying 
effective mass controls CL defined by the Green function. The static part + ■m'^/a{l — a) of 
the CL is connected with the longitudinal motion and is included in the Green function as well 
via the last phase shift factor (see Eq. (j^^ below). Consequently, the longitudinal momentum 
transfer is known and all the multiple interactions are included. 

The imaginary part of the LC potential Vqq{z2, r2, a) in Eq- (fTfj) is responsible for attenuation 
of the qq photon fluctuation in the medium, while the real part represents the interaction 
between the q and q. Because we are going to calculate the nuclear shadowing in DIS at 
medium and large one can safely neglect this q-q interaction as was done also in 

In the LC Green function approach jT3 CHI 113 CHI physical photon I7*) is decomposed 
into different Fock states, namely, the bare photon |7*)o, \qq), IqqG), etc. As we mentioned 
above the higher Fock states containing gluons describe the energy dependence of the photoab- 
sorption cross section on a nucleon. Besides, those Fock components lead also to GS as far 
as nuclear effects are concerned. However, these fluctuations are heavier and have a shorter 
coherence time (lifetime) than the lowest \qq) state. Therefore, at medium energies only \qq) 
fluctuations of the photon matter. Consequently, GS related to the higher Fock states will be 
dominated at high energies. Because we will calculate the nuclear shadowing at moderately 
small XBj (medium values of u) we can neglect the GS for our purposes. This is supported 
also by the main goal of this paper which is based on comparison of the realistic calculations 
for nuclear shadowing with the results obtained within the harmonic oscillatory Green function 
approach and presented in the paper [T3] where GS was neglected as well. 

One can describe a propagation of a noninteracting qq pair in a nuclear medium by the 
Green function satisfying the evolution Eq. (fTTjl . The LC potential Vqq{z2,r,a) in this case 
acquires only an imaginary part which represents absorption in the medium, 

ImVqq{z2, f, a) = p^{b, Z2) . (19) 

The analytical solution of Eq. ()17|1 is known only for the harmonic oscillator potential 
V{r-) oc r^. Consequently, one should use the dipole approximation 

aqq{r,s) = C{s)r\ (20) 

and uniform nuclear density 

PA{h.s)=p^Q{R\-h''-z'), (21) 

in order to to obtain the Green function in an analytical form. In Eq. fl21|) is the nuclear 
radius. The solution in this case is the harmonic oscillator Green function 



Gqq{r2,Z2]ri,Zi) = - — . exp<^ [r-^ + r^j cos{uj Az) - 2ri ■ r2 
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X 



exp 



ie^ Az 



a] 



where Az = Z2 — zi and 



where 



b{a) 



v a{l — a) ' 



(22) 



(23) 



h'^{a) = -ipA{h,z)va{l-a)C{s) . (24) 

The energy dependent factor C{s) in Eq. and the mean nuclear density po in Eq. (jTT]) 
can be adjusted by the procedure described in 16J. According this procedure the factor C(s) is 
fixed by demanding that calculations employing the approximation Eq. ()20j) reproduce correctly 
the results based on the realistic cross section (given by KST parametrization Eq. © or by 
GBW model Eq. (j^) in the high-energy limit Ic ^ Ra when the Green function takes the 
simple form (see Eq. (j3(Jj) below). Consequently, the factor C{s) is fixed by the relation, 



JdH Jd^r 


^5,''(r,a,Q2) 


2 ( 

1 1 — exp 


-lCis)r'T^ib) 


} 


J d'^r 


vl/J'^(f,a,g2) 


^C(s)r2 



JdH fd^r 


vl>^,^(f,a,Q2) 


2 r 

1 1 — exp 


-l(^qqir,s)TAib) 


} 


J d'^r 


^?i^(r,«,g^) 


2 



where 



TA{b) 



dz pA{b,z) 



(25) 



(26) 



is the nuclear thickness calculated with the realistic Wood-Saxon form of the nuclear density 
with parameters taken from [20] . This procedure is performed separately for T and L polarized 
photons and for each value of a. The value for the mean nuclear density po in Eq. (j2H) is 
determined in a similar way using relation 



dH 



1 - exp -cToPo jR\-b'^ 



dH 



1 — exp 



cToTAib) 



(27) 



The value of po turns out to be practically independent of the cross section ctq from 1 to 50 mb 
as was checked in jTHl 125] ■ 

We would like to emphasize that only in the high energy limit, Ic ^ Ra, it is possible 
to resum the whole multiple scattering series in an eikonal-formula. Correspondingly, the 
transverse separation r between q and q does not vary during propagation through the nucleus 
(Lorentz time dilation). Then the total photoabsorption cross section on a nucleus reads ^21\ 



a 



7*A, 



= 2 
= 2 



db Idr I da\'^gg{r,a,Q 



2\ |2 



exp 



■ (^qqir, s) TA{b) 



d'bll- 



exp 



-2 ^'?'?('^' s)TA{b) 



(2^ 



Note that the averaging of the whole exponential in (j^H|l makes this expression different from 
the Glauber eikonal approximation where a{r, s) is averaged in the exponent. 



Glauber i^^ Q 



2 d% 1 



exp 



7* AT, 



S,Q')TA{b) 



(29) 
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The difference is known as Gribov's inelastic corrections [HTj. In the case of DIS the Glauber 
approximation does not make sense (because of a small value of a'^''^, which is at most of the 
order of 100 fib for real photons) and the whole cross section is due to the inelastic shadowing. 

The eikonal formula ()28|) for the total photoabsorption cross section on a nucleus can be 
obtained as a limiting case of the Green function formalism. Indeed, in the high energy limit, 
z/ ^ oo the kinetic term in Eq. (fTTjl can be neglected and the Green function reads 



(5(r^ — fl) exp 



1 r^a 

Z J zi 



(30) 



After substitution of this expression into Eqs. (fT3j) and (fT^ one arrives at the result 

For smaller energies when ~ Ra, one has to take into account the variation of the 

transverse size r during propagation of the qq pair through the nucleus. This transverse size 

variation is naturally included using correct quantum-mechanical treatment based on the Green 

function formalism presented above. 

The overall total photoabsorption cross section on a nucleus is given as a sum of T and L 

polarizations, a'^*^ = aj^^ + e' a]^ ^, assuming that the photon polarization e' = 1. If one takes 

into account only qq Fock component of the photon the full expression after summation over 

all flavors, colors, helicities and spin states has the following form [221 



A I (fr I daagq{r,s){ '^^g{r,a,Q 



+ 



3 «p 



(27r)2 



Nf , 

^ r POO rco f f f 

Z'j Re / d^h / dzi / dz2 j da j d^ri / d^r2 

j^—l J J —oo J zi Jo J J 



X p^(6, zi) pAib, Z2) cTqq{r2, s) aqq{ri, s 

X 



e^!llIlK,{er,)K,{er2) 



+ 



ri r2 
2 „ 2 /1 „\2 



(31) 



mj + AQ'a'{l-ay K^^en) Ko{er2) Ggq{r-'2, Z2; n, z^] 



Here 



T,L 
qq I 



f, a, 



are the absolute squares of the LC wave functions for the qq fluctuation 



of T and L polarized photons summed over all flavors with the form given by Eqs. ()12|) and 
(fTSj) . respectively. 

In the high energy limit after substitution of expression (jHUj) for the Green function into 
Eq. (j31|) one arrives at the following results, which corresponds to Eq. (|28|) after inclusion of a 
sum of T and L polarizations : 



cr 



7* A 



{xBj,Q^ 



d% d\ 



If r 1 

(iaU-exp --aqq{r,s)TA{h) 



2 No a 



Nf 



C "em 



(2>r)- 



/=1 



a 



+ 



« + (1 - 



(32) 



a] 



Kl{er) 
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At photon polarization parameter e' = 1 the structure function ratio F^{xBj, Q"^) / {x Bj ■, Q"^) 
can be expressed via ratio of the total photoabsorption cross sections 



where the numerator on right-hand side (r.h.s.) is given by Eq. ()31|). whereas denominator can 
be expressed as the first term of Eq. divided by the mass number A. 

Finally we would like to emphasize that qq Fock component of the photon represents the 
higher twist shadowing correction [16j. This correction vanishes at large quark masses as 
l/rrtj. Not so for the higher Fock states containing gluons and leading to GS. GS represents 
the leading twist shadowing correction [331 • Besides, a steep energy dependence of the 
dipole cross section aqqir^s) (see Eqs. (0) and (jH))) especially at smaller dipole sizes r causes a 
steep energy rise of both corrections. 

3 Algorithm for numerical solution of the evolution equa- 
tion for the Green function 

As we mentioned in the previous section, an explicit analytical expression for the Green function 
Gqq{f*2^ 2:2; rl, can be found only for the quadratic form of the dipole cross section 

and for uniform nuclear density function (j2U)- It was already analyzed in JSJ that such 
an approximation should have a reasonable accuracy, especially for heavy nuclei. We also 
discussed in the previous Section that the higher accuracy can be achieved taking into account 
the fact that the expression (j28j) in the high energy limit can be easily calculated using realistic 
parametrizations of the dipole cross section (see Eqs. (0) and and a realistic nuclear density 
function p^(6, |30j. Consequently, one needs to know the full Green function only in the 
transition region from non-shadowing {xsj ~ 0.1) to a fully developed shadowing given by 
Eq. (j2Hl)- Therefore the value of the energy dependent factor C{s) in Eq. ((201) fixed [THl lTBj 
separately for T and L photon polarizations (see Eq. (|25p) in a such way that the asymptotic 
nuclear shadowing in DIS is the same for the realistic parametrizations of the dipole cross 
section Eqs. (0) and © and for approximation (pUj) . Correspondingly, the value po of the 
uniform nuclear density (PT|) was fixed in an analogical way as given by Eq. (j771) and described 
shortly in the previous Section. Such a procedure for determination of the factor C{s) and 
Po was applied also in \n\ ITH] with respect to incoherent and coherent production of vector 
mesons off nuclei. 

In order to remove the above mentioned uncertainties one should solve the evolution equation 
for the Green function numerically for arbitrary parametrization of the dipole cross section and 
for realistic nuclear density function. However, the tax for this general solution is that one does 
not obtain a nice analytical form for the Green function. First we present an algorithm for the 
exact numerical solution of the evolution equation. Using this algorithm we will calculate the 
nuclear shadowing in DIS and study how the new results change in comparison with predictions 
based on above mentioned approximations leading to harmonic oscillator Green function 



In the process of numerical solution of the Schrodinger equation (jl7j) for the Green function 
Gqq{f2i Z2;fl,zi) it Is uot Very convenient to treat the initial condition (jl8|) with two-dimensional 



Fi'{xBj,Q') 



4''ixBj,Q')+at''{xB„Q') 



(33) 



(El 
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Delta function on the r.h.s. In order to remove this problem one should use the following 
substitutions 



9i{r2, Z2; zi) = d ri Ko{eri) (Tgg(ri, s) Gqg{r2, Z2] Zi) 



and 



r2 , 
— 92\r2 
r2 



(fri Ki{eri) (7gg(ri, s) — Ggg{f*2, Z2; fl, zi 



(34) 



(35) 



Consequently, after some algebra with Eq. (fTTj) one can introduce new functions gi{r2, Z2; zi] 
and g2{r2, Z2', zi) which satisfy now the following evolution equations 

1 d ' 



. d 

dZ2 



2 IXqq 



dr2 r2 dr2 



and 



d 



1 



^-^92[r2, Z2; zi, 

dZ2 



92 1 d 



+ 



dri r2 dr2 



+ Vqq{z2,f2,a)\ gi{r2,Z2]Zi) (36) 



+ Vqq{z2, f2, ") \ 92{r2, Z2] Zi) , (37) 



[ 2 flqq 

with the boundary conditions 

gi{f2, Z2] zi%,=:,^ = Ko{er2)aqq{r2,s) (38) 

and 

5'2(^"2,^2;2i)|^2=zi = Ki{er2)aqq{r2,s) . (39) 
In Eqs. (jnni) and (jHTj) the quantity 

Hgq = ua {1 — a) (40) 

plays the role of the reduced mass of the qq pair. Consequently, the expression (|^ for total 
photoabsorption cross section on a nucleus now reads 



A I d^r I daaqq{r,s 



3 Op 



(2vr)^ 



Nf -, 

Zj Re / d^b / dzi / dz2 da d^r2 

J J —00 J Zl Jo J 



X pAib, zi) pA{h, Z2) aqq{r2, s) 



Ki{er2)g2{r2,Z2]Zi) 



(41) 



+ 



m] + AQ'^a^ (1 - af Ko{er2) gi{r2, Z2; zi] 



There are several approaches for solving the time-dependent one-dimensional Schrodinger 
equation (see for example jSH IS3 EEl EI] ) • One can not adopt directly these approaches for our 
purposes because one needs to treat the time-dependent two-dimensional Schrodinger equation 
(see Eq. fll7p ). Therefore, we will consider a modification of the method based on the Crank- 
Nicholson algorithm jHEj- Details of this method for numerical solution of Eqs. (j36|) and (j37jl 
are presented in the Appendix A. 
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4 Numerical results 



As we mentioned above the main goal of this paper is to present for the first time the reahstic 
predictions for nuclear shadowing based on exact numerical solutions of the evolution equation 
for the Green function. These predictions will be confronted with approximate results obtained 
using harmonic oscillatory form of the Green function (j221)- The proposed algorithm (see 
Appendix A) for numerical solution of the evolution equation for the Green function gives 
a possibility to calculate nuclear shadowing for arbitrary LC potential Vqq{z, r, a) and nuclear 
density function. It allows to perform an independent cross-check whether the results calculated 
using quadratic form of the imaginary part, ImVqq{z,r,a) = — 1/2 p^^b, z) C r"^ (C ~ 3), 
of the LC potential and uniform density function p^(6, 2) = po = 0.16 fm^^ correspond to 
predictions for nuclear shadowing taken from jl5j calculated using Green function of the form 
(j22I)- Therefore, we calculate nuclear shadowing for calcium and lead as was done in |15j . 

In order to realize one-to-one comparison with the results from the paper pTF] we made 
a several assumptions. As was mentioned in Section ^ we neglect the real part of the LC 
potential Vqq{z2,f2,C() in the Schrodinger equation (fTTjl ^ analyzing DIS at medium and large 
values of Q^. We neglect also the effects of nuclear antishadowing assuming that they are 
beyond the shadowing dynamics. The corresponding values of Bjorken xbj cover medium and 
medium large values, xbj G (0.001,0.1). For this reason we omit the effects of GS which 
are important at small xbj (large u). Although gluons can give some small (not negligible) 
contribution to nuclear shadowing at lower limit {xbj ~ 0.001) of investigated interval of xbj, 
for simplicity we do not include them in calculations as was done also in the paper ^21- 

We use an algorithm for numerical solution of the Schrodinger equation for the Green func- 
tion described in Section and Appendix A. Numerical solution of Schrodinger equation 
allows us to use realistic nuclear density function and realistic parametrizations of the dipole 
cross section, aqq{r,s). These parametrizations naturally incorporate the energy {xbj— ) de- 
pendence of aqq{r,s) which was not included so far in calculations (see [IS])^. We took the 
nuclear density function in Woods-Saxon form The realistic calculations of nuclear shad- 
owing were performed at two different parametrizations of the dipole cross section, GBW [21] 
given by Eq. ® and KST ^ given by Eq. ©. 

Nuclear shadowing effects were studied via xbj- behavior of the ratio of proton structure 
functions divided by the mass number A. The proton structure functions F2 {xBj,Q'^) 
and F-^i^XBjyQ"^) were calculated perturbatively (we fixed the quark masses at niq = 0.3 GeV, 
rris = 0.45 GeV and rric = 1.5 GeV) via total photoabsorption cross sections a^*^{xBj, Q"^) and 
(T'^*^(xBj, <5^) given by Eqs. (jH) and (jH)), respectively. The results of calculations are shown 
in Fig. 121 The dashed curves represent the predictions based on the harmonic oscillator Green 
function (j^^ approach corresponding to a constant nuclear density function po = 0.16 /m~^ 
and quadratic form of the dipole cross section (j^U)) with a constant factor C{s) ~ 3. The 
dotted curves correspond to calculations using the same quadratic form of the dipole cross 
section but realistic nuclear density function of the Woods-Saxon form jHU]. The thin and thick 
solid curve corresponds to realistic calculations based on the exact numerical solution of the 
evolution equation (fTTj) for the Green function using GBW (0) and KST (jUj) parametrization 

^Consequently, the real part of Vqq{z2, f2,a.) is neglected as well in differential equations and IpTTjl . 

^The energy dependence of dipole cross section was included only via energy-dependent factor C{s) in 
approximation H2U|) . The factor C (s) was determined by the procedure described shortly in Section [3 and 
presented in details in 1 161 1251 [T7] . 
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of the dipole cross section, respectively. 




0.8 - 



0.7 [ ■ — • — ■ — • — T . — . — ■ — . — 1 

0.001 0.01 0.1 0.01 0.1 

Figure 2: Nuclear shadowing for calcium and lead. The dashed curves 
are calculated using harmonic oscillator Green function approach (|^^ 
corresponding to constant nuclear density (j2H) and dipole cross section 
(j2(Jj) with C{s) ^3 [15 . The dotted curves are calculated for the same 
quadratic form of the dipole cross section but for realistic nuclear density 
function of the Woods-Saxon form 30J. The thin and thick solid curves 
correspond to exact numerical solution of the evolution equation for the 
Green function using GBW [23] and KST ^HI parametrization of the 
dipole cross section, respectively. 



At low XBj ^ 0.001 one should expect a saturation of nuclear shadowing at level given by 
Eqs. (j2Hl) or (jHH). For parametrization of the dipole cross section with constant C{s) ~ 3 
IjlSj , this saturation level is fixed at some value depending on and the nuclear mass number 
A (see the dashed and dotted lines in Fig. |21). However, it is not so for realistic parametrizations 
Eq. (jHI) and (jH} where the saturation level is not fixed exactly due to energy (Bjorken xbj-) 
dependence of the dipole cross section crqq{r, s). 
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In the process of realistic calculations of nuclear shadowing (|41|) we tested the correctness of 
tremendous computations based on the numerical evaluations of the functions gi and g2 from 
differential equations (jHUj) and ()37|) in a such way that in the high energy limit the results for 
nuclear shadowing must be the same as obtained from the expression p2|l . 

Results presented in Fig. El show quite a large deviation of the predictions within harmonic 
oscillator Green function approach (dashed lines) from realistic calculations performed for both 
parametrizations of the dipole cross section (thin and thick solid lines). This deviation de- 
pends on and the nuclear mass number A as a result of quadratic form of the dipole cross 
section (see Eq. (|20|) with C{s) ~ 3) and application of the constant nuclear density function, 
Po = 0.16 fm^'^. There is even not negligible difference between predictions using constant and 
realistic form of the nuclear density function (compare dashed and dotted lines). It allows to 
make conclusion that the form of nuclear density function is also important for model pre- 
dictions. Note, that the dashed curves correspond to predictions presented in the paper [THj. 
It is another cross-check for correctness of calculations using above presented algorithm for 
numerical solution of the evolution equation for the Green function. 

As one can see from Fig. |21 at small and medium values of {Q^ = 2 and 6GeV^), the 
approximate calculations depicted by the dashed lines agree better with realistic calculations 
using KST parametrization ^\ of the dipole cross section expressed by Eq. At large 

= 18 GeV^ and at xbj ^ 0.005, however, the dashed lines seem to be in better agreement 
with realistic calculations using GBW parametrization [25 given by Eq. (jj)). This fact confirms 
discussion presented in Section |21 that the GBW model is suited better at medium and large 

> 5 -h 10 GeV^ and at medium small and small xbj ^ 0.01 whereas the KST model prefers 
low and medium values of < 10 GeV^. 

Nevertheless, calculations of nuclear shadowing using harmonic oscillator Green function 
can be improved by determination of the energy- dependent factor C{s) in approximation ()20p 
by the procedure mentioned shortly above in Sect. |21 and described in details in fHl I2n| for 
calculation of nuclear shadowing in DIS and in f7| QHl for calculation of nuclear transparency 
for coherent and incoherent vector meson production off nuclei. That procedure allows to 
evaluate the factor C{s) for each cm. energy squared s depending on the values of and A. 
As a result, at fixed the parameter C{s) rises with s as a consequence of energy- {xbj- ) 
dependent realistic dipole cross section given by Eqs. (0) and Thus, the value of C{s) at s 
corresponding to xbj G (0.001,0.01) exceeds so the fixed value C{s) ~ 3 used in predictions in 
[15, . This fact should lead to a larger nuclear shadowing at xbj G (0.001,0.01) in comparison 
with what is shown in Fig. |21 by the dashed lines. 

As we mentioned above in Sectional the difference between the thin and thick solid lines in 
Fig. 121 can be treated as a measure of the theoretical uncertainty in the kinematic region where 
the both realistic parametrizations GBW (0) and KST © are applicable. Therefore, it would be 
very useful for the future realistic calculations to connect advantages of both parametrizations 
in the modified model for dipole cross section which can be then safely used for all dipole sizes 
covering perturbative as well as nonperturbative region. 

5 Summary and conclusions 

We present a rigorous quantum-mechanical approach based on the light-cone QCD Green func- 
tion formalism which naturally incorporates the interference effects of CT and CL. Within this 
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approach [1^1 1201 HE] we study nuclear shadowing in deep-inelastic scattering at moderately 
small Bjorken Xbj- 

Calculations of nuclear shadowing performed so far were based only on the efforts to solve 
the evolution equation for the Green function analytically. Analytical harmonic oscillatory form 
of the Green function could be obtained only taking into account additional approximations 
like a constant nuclear density function 1)211) and the dipole cross section of the quadratic form 
(j2ni)- It brings additional theoretical uncertainties in predictions for nuclear shadowing. In 
order to remove these uncertainties we solve the evolution equation for the Green function 
numerically. 

We perform for the first time the exact numerical solution of the evolution equation for 
the Green function using two realistic parametrizations of the dipole cross section (GBW |21] 
and KST [^) and realistic nuclear density function of the Woods-Saxon form [30]. This exact 
numerical solution does not require to put any additional approximations. Analyzing only 
medium and large values of we neglect the real part of the LC potential Vqq{z2, r^, a) in the 
time-dependent two-dimensional Schrodinger equation (|T7|l responsible for interaction between 
q and q. We neglect also the nuclear antishadowing effect as was done in [12] assuming that it is 
beyond the shadowing dynamics. Performing calculations at medium and medium large values 
of Xsj G (0.001,0.1) we neglect for simplicity also the contribution of the higher Fock states 
leading to effects of GS. This is supported also by the one-to-one comparison of the realistic 
calculations with the predictions from the paper (THj, where GS is neglected as well. 

In order to compare the realistic calculations with data on nuclear shadowing, the effects of 
GS should be taken into account especially at xbj ^ 0.001. The same path integral technique 
[in] can be applied in this case. However, the calculations of GS (see jSH], for example) were 
performed so far using analogical approximations as already mentioned above like a constant 
nuclear density function and the quadratic form (see Eq. ((201)) of the dipole gluon-gluon-nucleon 
cross section, Cgglr, s) = 9/4agq{r, s) plus further assumptions which simplify the final expres- 
sion for GS. Moreover, the GS was calculated from the shadowing of the \qqG) Fock component 
of a longitudinally polarized photon at sufficiently large where the three-body Green func- 
tion, GqqG, is assumed to be factorized as a product of two-body ones Using the algorithm 
presented above one can calculate GS exactly for the general case of nuclear shadowing for a 
three-parton system, i.e. one can solve numerically the Schrodinger equation for the Green 
function Ggqc describing propagation of the qqC system through a nuclear medium. We are 
going to calculate numerically the gluon contribution to nuclear shadowing in a forthcoming 
paper. 

We present analogical numerical results of nuclear shadowing in DIS with correct quantum- 
mechanical treatment of multiple interaction of the virtual photon fluctuations and of the 
nuclear form factor as was done in We found quite large differences (see Fig. ^ 

between realistic predictions and the approximate results obtained within harmonic oscillator 
Green function approach (j221)- At small and medium values of the approximate predictions 
agree better with realistic calculations using KST parametrization ^H] of the dipole cross section 
Eq. (jni). At large Q^, however, they seem to be in better agreement with realistic calculations 
using GBW parametrization [25 given by Eq. 1^. It confirms the fact that the GBW model 
is well suited at medium and large and at medium small and small xbj, whereas the KST 
model prefers low and medium values of Q^. Therefore, the future realistic calculations require 
to revise existing parametrizations for dipole cross section in order to be used for whole region 
of dipole sizes. 
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Concluding, the universality of the LC dipole approach based on the Green function formal- 
ism allows us to apply the presented algorithm for the exact numerical solution of the evolution 
equation for the Green function also for calculations of other processes like Drell-Yan produc- 
tion, vector meson production etc. including the effects of gluon shadowing at high energies as 
well. 

Acknowledgments: We are grateful to Alexander Tarasov for stimulating discussions. This 
work has been supported in part by the Slovak Funding Agency, Grant No. 2/2099/22 and 
Grant No. 2/1169/21. 

Appendix A Description of the method for numerical 

solution of the time-dependent Schrodinger 
equation 

We treat here only differential equation pHjl for the function gi{f2, Z2]Zi) and describe in details 
the method for its numerical solution. This method is then analogically applicable for numerical 
solution of Eq. (|n7j) for the function (72(^^5 -22; -^i)- 

Looking at Eq. (IHUj). one needs to solve numerically the following time-dependent Schrodinger 
equation ^ 

ij^g^{f,t)=Hg^{f,t), (A.l) 

where H is the Hamiltonian operator defined by 

1 



H 



2 l^qq 



+ Vqq{r,t). (A.2) 



Here the complex LC potential Vqq{r, t) is assumed to have only the imaginary part responsible 
for absorption of qq photon fluctuation in the nuclear medium (see discussion in Section |21 and 
Eq. (HHD), 

Vqqir, = (^qg{r, s) pA^b, t) , (A.3) 

where b is the nuclear impact parameter and r is the transverse separation between q and q at 
the point Z2- The longitudinal coordinate Z2 plays the role of time t for the qq pair propagation 
from the point zi. In Eq. (1A.2|1 the quantity is expressed by Eq. pn|l . 

Ignoring for a moment the fact that H is an operator, Eq. (jA.ljl has the formal solution, 

gi{f,t) = exp{-ttH)gi{f,0), (A.4) 

where 5^1 (r, 0) is the function at t = 0. Thus, if one knows gi{f, 0), one can formally calculate the 
behavior at all future times using Eq. ()A.4|1 . Unfortunately, this formal solution is not of much 
practical use since the Taylor expansion of the exponential factor in Eq. ()A.4j) involves a very 
large (infinite) number of terms. However, it does suggest a way to proceed numerically. Let 
us consider a formal solution applying over a very small time interval. After time discretization 
in steps At using Eq. ()A.4|) we obtain 

gi{f,t + At) = exp{-iAtH)gi{f,t). (A.5) 



^We put r2 = r and Z2 plays the role of time t. 
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Consequently, at sufficiently small time intervals At the higher order terms in Taylor expansion 
of the exponential factor in ()A.4|) are small enough and can be neglected. Then we can include 
only the term linear in 

exp{-iAtH)^l-iAtH. (A.6) 

However, this way of approximating the exponential factor is not correct with regard to main- 
taining unitarity [SHI- order to settle this problem one should use an approach that satisfy 
unitarity writing the exponential factor in Eq. ()A.5|1 in what is known as the Cayley form 



. , 1-UAtH 

exp (-i At H) ^ ^ ^ 

1 + UAtH 



(A.7) 



Using this approximation and Eq. ()A.4|1 to propagate the function gi{f,t) forward in time we 
obtain 

1-UAtH , ^ , , 

g^{f, t + At)^ f-— ^ (71 (f, t) . (A.8) 



1 + liAtH 



This expression will be the basis for numerical approach. From ()A.8jl we ffist obtain 



l + -iAtH 
2 



g^{f,t + At) 



1- - iAtH 
2 



giir,t) . 



(A.9) 



Given flA.9|) . a natural way to proceed is to discretize also space into units of size Ar and 
write the function gi as gi{r,t) = gi{mAr,nAt). Then one can express the ffist and second 
derivatives included in the Hamiltonian operator ()A.2|) in the usual finite-difference form : 



d gi{r + Ar) - gi{r,t) gi{m + l,n) - gi{m,n) 

— gi{r,t) ^ = 

o r Ar Ar 



(A.IO) 



and 



dr' 



■9i{r,t) 



gi{r + Ar) - gi{r,t) 
Ar 



d_ 

dr 

gi{m+ l,n) -2gi{m,n) + gi{m - 



gi{r + Ar,t)~2gi{r,t)+gi{r-Ar,t) 



(Ar) 



(Ar) 



(A.ll) 



If one replaces in Eq. ()A.9|) the Hamiltonian operator by ()A.2|) . converting everything to 
finite-deference form using also Eqs. ()A.10|) and ()A.11|) . and rearranging a few terms one obtains 
the following expression 



(7i(r/i + 1, n + 1) + h{m) 



2i\-2figg {Arf V{m, n + 1) - ( Ar) 



1 

2 

m 



gi{m,n + 1) 



+ h{m) gi{m — l,n + 1) 



= —gi{m + l,n) + h(m) 

— h{m) gi{m — l,n) 
where the function h{m) is given by 



2iX + 2figg (Ar)2 V{m, n) + {Arf + 2 + 



gi[m,n) 



h{m) 



m 

1 + m 



(A.12) 



(A.13) 
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and 



A 



2(Ar)V,-, 
At 



(A.14) 



with the reduced mass of qq pair defined by Eq. (pH) . 

The algorithm very effective for solving the time-dependent Schrodinger equation in one 
dimension is known as the Crank-Nicholson method described in details in [3111101 ■ However, 
in order to solve ()A.12|) one should modify this method for more complicated case of two 
dimensional Schrodinger equation. We begin by defining a shorthand for the r.h.s. of Eq. ()A.12|) 



Q['m,n) = —g 1(171 + l,n) + him) 

— /i(m) (7i(m — 1, , 
in order to rewrite Eq. ()A.12|) as 



2 z A + 2 /iq-g {Arf V{m, n) + {Arf + 2 + 



m 



gi[m,n) 
(A.15) 



gi{m + l,n + l) + h{m) 



2i\-2fxgg (Arf V{m, n + 1) - {Arf - 2 



m 



+ /i(m) 5fi(m — 1, n + 1) = Q{m,n) . 



gi{m,n + 1) 
(A.16) 



For convenient numerical procedure one should write gi{m + l,n + 1) as a function of just 
(7i(m, n + 1) in the following form 



gi{m + l,n + 1) = e(m, n) gi{m, n + 1) + /(m, n) 



(A.17) 



and so one can calculate gi{m + l,n + 1) directly from gi{m, n + 1). 

If one inserts Eq. flA.17|) into Eq. ()A.16|) and does a little arithmetic, one can find that the 
factors e(m, n) and f{m,n) must be given by the following implicit relations 



e m, n) 



h{m) 



2 + - + 2figg {Arf V{m, n + 1) + {Arf - 2 i A - — r 

m eim — l,n) 



and 



/(m, n) = f2(m, n) + /i(m) 



/(m - l,n) 
e(m — 1, n) 



(A.18) 



(A.19) 



One supposes that a propagation of qq pair in the nuclear medium is confined to some region 
of space so that spatial index runs from m = to m = M and imposes the boundary conditions 
gi{0,n) = gi{M,n) = 0. The expressions ()A.18|) and ()A.19|) for the factors e{m,n) and f{m,n) 
can be applied only in the interior of the system. Consequently, from the boundary condition 
for the function gi at m = together with Eqs. ()A.15|1 and ()A.17|1 one can find that at this end 
of the system 



e(l,n) = h{l) 



and 



3 + 2 fxgg {ArY V{l,n + 1)+ (Ar) 



f{l,n) = nil,n). 



2iX 



(A.20) 



(A.21) 



For the first time step, n = 0, the factors fi(l,0), e(l,0) and /(1,0) can be explicitly 
calculated from the initial function, gi{m,0) = aqq{m Ar, s) Ko^em Ar), which is assumed to 
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be given as an initial condition (see also Eq. (jHHj) ). Using known values of e(l,0) and /(1,0) 
one can calculate then the factors e(2, 0) and /(2,0) from the implicit expressions Eqs. ()A.18|) 
and ()A.19|) and continue so for all values of m along the system. Hence, we traverse the system 
from m = to m = M, to calculate e(m, 0) and f{m, 0) for all m. 

For the further purposes, Eq. ()A.17|) can be rearranged in the following form 

/ , gi{m + l,n+l) - /(m,n) 

gi{m,n + l) = . (A.22) 

e(m, n) 

As was mentioned above, at the end of the system m = M the function gi vanishes. Conse- 
quently, one can write 

(M 1 9iiM,n + l)-fiM-l,n) /(M-l,n) 

gi{M -l,n+l) = — — - — = -__ — - — - (A.23) 

e(M — l,n) e(M— l,n) 

since gi{M,n) = for all values of n. One can thus use Eq. ()A.23|) to obtain gi{M — 1, 1), 
which is the value of the new function one spatial unit in from the "right" boundary, m = M. 
Then one can use Eq. flA.22|) to calculate the function gi at m = M — 2, M — 3, etc., as one 
traverses the system backward from large to small values of m. 

Finally, the algorithm can be summarized in the following way (see also the schematic 
description in Fig. |2)) : 

1. One begins with the initial function gi{m, 0) given by Eq. (|HHj). 

2. The system is traversed from small to large values of m and the functions e(m, 0) and 
f{m, 0) are calculated using Eqs. (IA.20jl and (IA.21I1 initially and Eqs. flA.lSjl and (IA.19j) 
thereafter. 

3. The system is traversed from large to small values of m and gi{m, 1) is calculated using 
Eq. ()A.23|) initially and Eq. ()A.22|) thereafter. This completes one iteration and yields 
the function at n = 1 (t = At). 

4. Steps (2) and (3) are repeated to obtain the function gi as a function of time {n > 1). 
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